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We study the lowest order conservation laws in one-dimensional (ID) integrable quantum many- 
body models (IQM) as the Heisenberg spin 1/2 chain, the Hubbard and t-J model. We show that 
the energy current is closely related to the first conservation law in these models and therefore 
the thermal transport coefficients are anomalous. Using an inequality on the time decay of current 
correlations we show how the existence of conserved quantities implies a finite charge stiffness (weight 
of the zero frequency component of the conductivity) and so ideal conductivity at finite temperatures. 
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Finite temperature transport properties, as electrical Q] and thermal Q conductivity or spin dynamics j^] of 
strongly correlated ID are recently in the center of experimental and theoretical interest. In this context, IQM's as 
the Heisenberg spin 1/2 chain, the Hubbard or t — J model B are often the starting point for describing ID systems. 
Further, from another point of view, the integrability of these models offers the possibility of studying the interplay 
between chaotic properties and transport in quantum many-body systems, mostly studied in the context of classical 
many-body Q or single particle quantum problems . 

Based on analytical and numerical studies, we have recently proposed that IQM show anomalous, dissipationlcss 
finite temperature conductivity [^-[n^. It is natural to think that this anomalous transport behavior is related to 
the macroscopic number of conserved quantities JIT|-p^| characterizing these systems. A set of conservation laws is 
represented by local involutive operators Q n , commuting with each other [Q n ,Q m ] = and with the Hamiltonian, 
[Qni H] = 0. The index n indicates that the operator Q n is of the form Q n = where are local operators 

involving n sites around site i, on a lattice of L sites. 



Although rather formal procedures exist for the construction of these operators [g_3 14 1, it is not clear how to study 
their physical content and even more how to take them into account in the analysis of transport properties. In this 
paper we show, for different models of actual interest, that already the first nontrivial quantity Q3 (Q2 often denotes 
the Hamiltonian) has a simple physical significance: it is (or it is closely related to) the energy current operator. 
Further, we analyze how the coupling of the energy current or current operator to the conserved quantities results to 
time correlations not decaying to zero at long times. Thus, transport has not simple diffusive character and within 
the Kubo linear response theory [fl5f is described by diverging or ill-defined transport coefficients. 

We can relate the time decay of correlations to the local conserved quantities in the Hamiltonian systems we discuss 
by using an inequality proposed by Mazur Jig] : 

lim I f T < A (t)A >dt>Y < AQ " >2 , (1) 

where <> denotes thermodynamic average, the sum is over a subset of conserved quantities Q n , orthogonal to each 
other < QnQm >=< Q n > <5n,m, = A and < A >= 0. In the following, we will only consider Q3 in relation (|l|) 
so that the issue of orthogonality will not enter. Further, we can write < A(t)A >= Caa + C(t) as the sum of a 
time-independent factor, 

Caa = Y.P- E l<«l^>| 2 (2) 

a b(c b =e a ) 



and a time-dependent one, 



C(t) = J2Pa E < a\A\b > IV^"^*. (3) 



Here \a >, \b > arc cigcnstates of the Hamiltonian, p a = e /3c ° /Z the corresponding Boltzmann weights and f3 the 
inverse of the temperature. For time correlations < A(t)A > with non-singular low frequency behavior, the term 

liniT^oo y J a G(t)dt goes to zero and so Caa = lim^oo < A(t)A >, 
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In particular, we will use this inequality in the analysis of the real part of the conductivity, a'(uo) — 2ttD(T)S(uj) + 
<j reg (u>), related, within linear response theory, to the current-current correlation < J(t)J >. A finite value of the 
charge stiffness £>, given also by D = \<jJ<t" '(w)| w _>0j implies an ideally conducting system [|l7|,||. We will now argue 
that D ~ 2L^JJ an d therefore the following inequality holds for the charge stiffness: 



In this derivation, we assume again that the regular part of the conductivity a reg {io) shows a non-singular behavior at 
low frequencies so that the contribution from C(t) i n (|l| ) vanishes. This is a very mild condition for the physical systems 
we consider; actually, we even have indications P,|10(| from numerical simulations, that these IQM are characterized 
by a pseudogap and so a vanishing regular part o~ reg (u> — > 0) . 

To relate D to Cjj, it is convenient to use a generalization of Kohn's approach |l7]J^] to finite temperatures. In this 
formulation D(T) can be expressed as the thermal average of curvatures of energy levels in a Hamiltonian describing 
a system on a ring pierced by a fictitious flux (f>, D = J2 a Pa ( iJ^f" ) U-»o ■ Evaluating the second derivative of the 
free energy Fasa function of the flux <fi we find: 

a a T 

In the models we will discuss, the third term on the r.h.s. vanishes by symmetry (summing over k and —k subspaces). 
Further, these systems show no persistent currents at finite temperatures in the thermodynamic limit, therefore 
T^p-U— >o — > in this limit. We have numerically verified that this is indeed the case for temperatures larger than 
the level spacing; at zero temperature, there is no contradiction with Kohn's formula for D as the second term in the 
r.h.s. vanishes in a ground state with zero current. Finally, as ^l^o =< °lil a > an d degenerate levels contribute 
a vanishing weight in (|J), we arrive at relation ([5]). This inequality provides the lower bound for the charge stiffness 
D which, if not zero, implies ideal conductivity at finite temperatures. 

In general, it is difficult to evaluate the right hand side of the inequality (|]) involving the "overlap" < AQ n >. 
However, we will give some examples below, this correlation can easily be evaluated in the case of a grand canonical 
trace over states, in the thermodynamic limit and for j3 — ► 0. We thus obtain the charge stiffness in leading order in 
(3. 

Before studying concrete models, we construct the energy current operator as follows: we consider Hamiltonians 
defined on a ID lattice with L sites and periodic boundary conditions /il,z,+i = of the form: 

L 

= (7) 

i=l 

Because the energy is a conserved quantity, the time evolution of the local energy operator hi^\{t) can be written 
as the discrete divergence of the energy current operator J E — Yli=i if'- 

dhi t i + i{t) . rrT L ^ E , E 



dt 



= i[H,h i>i+1 (t)] = -(jt+i(t)-3?(t)), (8) 



where h iii+ i{t) = e' lHt h^ i+ ie lHt and jf = — Now, by direct evaluation for several IQM, we will sec 
that the energy current J E is closely related to a conserved quantity. 

i) Heisenberg model: The general anisotropic Heisenberg Hamiltonian is given by: 

L L 

H = = yXJgfff Sf +l + J y SfS^ +1 + J z Sf (9) 

1=1 t=l 

where Sf — \af , af are Pauli spin operators with components a — x, y, z at site i. The local energy current operator 
if is: 

jf = J x J y {xzy - yzx)i- lti+ i + J y J z {yxz - zxy)i-i,i+i + J z J x {zyx - xyz)i-i, i+ i, (10) 
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with (a/J 7 - 7j8a)i-i,i+i = Sf-i^+i " S?-iS?S& v 

Now it is straigthforward to verify that the global energy current operator, J £ = ^2f =1 jf, commutes with the 
Hamiltonian (j^) and coincides with the first nontrivial conserved quantity, usually denoted by Q3, obtained from an 
expansion of the transfer matrix in the algebraic Bethe ansatz method jll],|l4| . In agreement with the notation Q2 for 
the Hamiltonian, the local energy operator hi i+i involves two sites + 1), while the local energy current operator 
qf = j E involves three sites (i — 1, i, i + 1). 

The vanishing commutator [J , H ] = implies that the energy current time correlations are independent of time. 

< J E (t)J E >=5> q jL, (11) 

a 

where js a are the eigenvalues of J E , J E \a >= jEa\a >• The non-decaying of the energy current leads to a diverging 
thermal conductivity related to the < J E (t)J E > correlation flq| . 

As for the conductivity, it is more relevant to discuss the fermionic version of the Heisenberg model, defined through 
the Jordan- Wigner transformation, the so called t — V model: 

L L 

H = (-t)J2(4^+i + h.c) + Vj2( n * ~ ^K+i - 2 } ' (12) 

i=l i—1 

where Cj(cj) denote annihilation (creation) operators of spinless fermions at site i and = c\c{. 

In this case, the corresponding energy current operator that commutes with the Hamiltonian ( |l2| ) is given by: 

J E = ^(-t) 2 (ict +1 Ci_i + h.c.) + Vj^i+xirn-x + n l+2 - 1), (13) 

i 

where ji.i+i = ( — ^)(~* c l+i c i + h.c.) is the particle current. Therefore, for this fermionic model, we find that the 
< J E (t)J E > as well as the < J E (t)J > correlations are time independent implying a diverging thermal conductivity 
and ill-defined thermopower respectively. 

Regarding the charge stiffness D, we can evaluate analytically < JQ3 > 2 / < Q\ > for — > and in the 
thermodynamic limit, obtaining from ([5]), 

D> P_ 2V* P {l- P ){2 P -iy 

- 2 l + V 2 {2p 2 -2p+l)' [ ' 

where p is the fermion density. We note that for p ^ 1/2, D is finite implying ideal conductivity as we have suggested 
before |^]. For p = 1/2, this inequality is however insufficient for proving that D is nonzero. Due to the electron - hole 
symmetry, this remains true even if we consider all the higher order conserved quantities Q„ - The reason is that, for 
the Heisenberg model, all Q n 's can be generated [||[|4) by a recursive relation: [B, Q n ] ~ Q n +i where B is a "boost" 
operator given by: B — ^2 n nh ntn -\-%. Then by the electron - hole transformation Cj = ( — 1 ) 8 c| , we see that J — > — J 
but Q n — > Q n and therefore for p = 1/2, < JQ n >= 0. The eventual nonorthogonality of the Q' n s is not important 
as we can see by considering new orthogonal conserved quantities constructed using, for instance, a Gram-Schmidt 
orthogonalization procedure. 

In Table |, we present some indicative numerical results comparing Cjj with < JQ3 > 2 / < Q§ > for a couple of 
p values and j3 — > 0. The results for Cjj were obtained by exact diagonalization of the Hamiltonian matrix on finite 
size lattices [L up to 20 sites), followed by finite size scaling using a 2nd order polynomial in 1/L. From this table 
we see that: i) the smaller the density the more the inequality @ is exhausted by just considering the contribution 
from Q3, ii) for V/t — » 00 the overlap < JQ3 > gives the total weight of Cjj; indeed, studying the higher order local 
conserved quantities we can see that they only contribute terms in powers of 1/V. Nevertheless, it is not clear why 
the inequality (Q) is exhausted and no other, e.g. nonlocal, conserved quantities contribute. 

Finally, returning to the Heisenberg model, we note that the bound (|l4|) implies anomalous spin transport at finite 
magnetization, as the particle current maps to the spin current and the density p to the magnetization (half-filling 
p = 1/2 corresponding to zero magnetization). 

ii) Hubbard model: It describes a system of interacting fermions on a lattice with Hamiltonian given by: 

L L 

H = (-*) ]T (c.U +lff + h.c.) + t/J> 4T - -)(n iL - -), (15) 
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where Ci CT (c| CT ) are annihilation (creation) operators of fermions with spin a =|, | at site i and — c\ a Ci a . 
Similarly as above, we can define a local energy operator by: 

Ki+i = M)^( c ^ c *+ lff + h - c -) + 2"(K't - \)i n ii ~\) + - ~ \))- ( 16 ) 

From the time evolution of /ij^+i we find the local energy current operator jf involving sites [i — 1, i, i + 1): 

if = ^(-tfiicl+laCi-lv + h - C -) - ^-(ji-hhV + ii,i+l,cr)K,-<7 - (17) 



where ji^+ia = (—t)(—icl + -,„Ci C + h.c.) is the particle current. By comparing this expression for the energy current 
to the conserved quantity p2|Jl4| Q 3 , we find that they coincide when the factor U/2 in ( |l7| ) is replaced by U. So 
the energy current J E — X)i=i if does not commute with the Hamiltonian. However, as J E has a finite overlap 
< J E Qz >, with the conserved quantity Q3 we still find that the energy current correlations decay to a finite value 
at long times so that the thermal transport coefficients are anomalous. We can find a lower bound for the decay by 
using (uh for (3 — > and in the thermodynamic limit: 



< J E Q 3 > 2 
A~ — '--■>-•>- '- < Qi > 



lim < J E (t)J E >= Cjeje > (18) 



<J B Q 3 > 2 _ > , U 4 E g 2p g (l-p g )(2,£ ff -2p-<7+l)] 2 

<QI> + 4 £ CT 2p a (l-p CT )[l + [/W- ff -2p-. + l)r (19) 

As for the charge stiffness D, we can again evaluate analytically < JQ3 > 2 / < Q§ > for (3 — > and in the 
thermodynamic limit, obtaining from ([5]), 

D > g [£/E.2p g (l-p CT )(2p_ CT -l)] 2 

" 2 E CT 2^(1 - P<r )[l + C/ 2 (2p 2 _ c7 - 2p_ ff + 1)] ' 1 J 

where p a are the densities of a =T, J, fermions. For p a = 1/2, the right hand side of ( p0| ) vanishes, although a general 
proof involving all higher conserved quantities is not possible as a boost operator for the Hubbard model is not known. 

Hi) "t-J" model: It belongs to a class of multicomponent quantum systems |fl8f describing interacting particles of 
different species, singly occupying each site. The Hamiltonian acts on each bond + by the operator P-i^+i which 
permutes neighboring particles, independently of their type: 



H = J2^+i (21) 

»=1 

For this generic model we can directly verify that the energy current operator 

L 

J B = -i^2[Pi-i,i,Pi,i+i], (22) 
i=i 

coincides with a conserved quantity fl9| and so commutes with the Hamiltonian. 

Now considering three types of particles, corresponding to empty sites, up spins and down spins, we recover the 
t — J model (l9) for special values of J/t. This model describes a system of interacting fermions subject to a constraint 
of no double occupancy, with Hamiltonian given by: 

L L 

H = -tJ2 p (4*ci+i° + h.c)P + Jj20i&+i ~ «<n<+i/4) +2JV - i, (23) 

(7,2 — 1 i — 1 

where c,- CT (ct) are annihilation (creation) operators of a fermion on site i with spin a =|, {. P = n^=i(l — n i~\ n il) 1S 
a projection operator on sites with no double occupancy, rii a — c\ a Ci ai N — J^ i=1 a ni a . 
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This model is integrable for J ft = 0, corresponding to the U/t — > oo limit of the Hubbard model or to the model ( |2~l| ) 
where permutations act only on bonds with "empty - up spin" or "empty - down" configurations. For this case, we 
found that the corresponding energy current commutes with the Hamiltonian, as is also known for the particle current 
pof . Finally, for J = 2t, the "supersymmetric" model ( p3| ) is also integrable and the energy current coincides with the 
conserved quantity Q3 as is presented in reference (19). Therefore the transport coefficients of the supersymmetric 
t — J model are also anomalous. 



The above results imply that, at least, certain quantities related to transport coefficients in IQM are non-ergodic 
(see however reference H for the recent notion of "mixing"). Within linear response theory, this translates to ideal 
conducting behavior at finite temperatures, the charge stiffness D being a measure of nonergodicity. We also expect 
that noise, as is described by the spectral properties of the current-current correlations, shows anomalous behavior, 
more characteristic of a ballistic rather than a diffusive system. This behavior is to be contrasted to the normal 
dissipative behavior we found ]9|,[To| for similar nonintegrable systems with no conservation laws. So IQM, on the 
one hand are not generic models for studying finite temperature transport in quantum many-body systems, but on 
the other hand they offer the possibility of observing new effects. The main remaining question of course, is the 
robustness of (nearly) ideal conducting behavior for systems close to integrability, a problem very similar to the one 
in classical near-integrable nonlinear systems. Anyway, as it is known that conservation laws play an important role 
on transport properties [^lf they should be taken into account in approximate or exact analysis of these properties. 
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TABLE I. < ®% > /°JJ as a function of filling p and V/t. 
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